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Abstract. We find the transition kernels for four Markovian interacting particle systems 
' on the line, by proving that each of these kernels is intertwined with a Karlin-McGregor 

, type kernel. The resulting kernels all inherit the determinantal structure from the Karlin- 

' McGregor formula, and have a similar form to Schiitz's kernel for the totally asymmetric 

04 ' simple exclusion process. 

o 



Resume. Nous trouvons les noyaux de transition de quatre systemes markoviens de partic- 
ules en interaction sur une ligne, en prouvant que chacun de ces noyaux s'entrelace avec un 
noyau du type de Karlin-McGregor. Tons les noyaux resultants heritent de la structure de 
determinant de la formule de Karlin-McGregor et ont une forme similaire acelle du noyau 
de Schiitz pour le processus d'exclusion simple totalement asymetrique. 



r-| '. 1. Introduction 

» . 

j3 ' Non-colliding Markov processes are canonical examples of stochastic processes with a 

determinantal transition kernel, given by the Karlin-McGregor formula. A determinantal 
transition kernel of a different form, yet similar to the Karlin-McGregor kernel, was encoun- 
CN ■ tered by Schiitz [IT] in his study of the totally asymmetric simple exclusion process. This 

^ ! work has stimulated much recent research, e.g., [11 \T0 \ [T ^ \W [ [2T | [22]. 

I In this note we explicitly connect Schiitz type formulae for particle systems and Karlin- 

OO ■ McGregor type formulae for non-colliding processes. Our approach builds upon deep con- 

nections between particle processes and non-colliding processes (or random-matrix theory), 
' which have been recently discovered [U |6l [9l [13] . A combinatorial correspondence known as 

, the Robinson-Schensted-Knuth (RSK) correspondence links these processes. This RSK cor- 

^ I respondence gives a coupling of the non-colliding process and the particle process. Our key 

contributions are that this coupling implies an intertwining of their transition semigroups, 
and that this intertwining is enough to find the transition kernel of the particle process. 
^ ' We give a number of new formulae by applying this method to four variants of the RSK 

^ . correspondence. Augmented with systems arising from suitable limiting procedures, the 

particle systems we treat in this way are known to play pivotal roles in a wide range of 
interesting applied problems. For instance, they appear in the context of queues in series, 
last-passage percolation, growth models, and fragmentation models. 

This note is organized as follows. In Section [21 we introduce four interacting particle 
systems and we present the associated transition kernels. Section [3] describes the four 
variants of the RSK correspondence we use in our analysis, and derives the aforementioned 
intertwining of the semigroups. Finally, it is the topic Section [D to use this intertwining for 
finding the transition mechanism of the interacting particles. 

2. Interacting particles on the line; main results 

We are concerned with a system of N particles, each with a position in the integer lattice 
Z, evolving in discrete time. We will consider four possible cases. Particles will move from 
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the left to the right making either Bernoulh or geometrically sized jumps, with one of two 
possible interactions that maintains their relative orderings (blocking or pushing). We begin 
by describing these four processes more precisely. In each case Yi{n) denotes the position of 
particle number i at time n. We order the particles, so that Y takes values in either W'^ or 
, where 

= {z eZ^;zN < ZN-i < ... < zi) 
= {z eZ^;zi < Z2 < ... < zn). 
Throughout, we let p = {pi,P2, ■ ■ ■ ,Pn) be a vector with each pk G (0, 1). 

CASE A: Geometric jumps with pushing. Particles are labelled from left to right, 
so Yi{n) < Y2{n) < . . . < YN{n). Between time n — 1 and n, each of the particles moves to 
the right according to some geometrically distributed jump, having parameter pi for particle 
i. The order in which the particles jump is given by their labels, so the leftmost particle 
jumps first. Overtaken particles (if any) are moved to the same position as the jumping 
particle, a position from which the next particle subsequently makes its own jump. One can 
thus think of particles 'pushing' other particles to maintain their relative orderings. 

This leads to the following stochastic recursion. The evolution is generated from a family 
(^(fc,n);A: G {1, 2, . . . A^}, n G N) of independent geometric random variables satisfying 
P{^{k, n) = r) = (1 —pkjPk for r = 0, 1, 2, . . ., via the recursions Yi{n) = Yi(n — 1) +C(1, n), 
and for A; = 2, 3, . . . N, 

Yk{n) = max{Yk{n- l),Yk^i{n)) +^{k,n). 

Note that Y = {Y{n); n > O) is a Markov chain on . 

One application area where this recursion arises is the theory of queueing networks. 
Indeed, the particle system with exponentially distributed jumps, which is obtained after a 
suitable limiting procedure, relates to a series Jackson network. Here Y{n) corresponds to 
the departure instants of the rzth customer from each of queues in series. These networks 
are investigated further in our companion paper [5]. 

The vector Y also plays an important role in the context of directed last-passage per- 
colation with geometrically distributed travel times and 'origin' (1,1), where Y{n) can be 
interpreted as the vector of maximal travel times to the sites (n + 1, 1), . . . , (n -f- 1, N). Very 
recently, Johansson [10] has derived the transition kernel of Y in the case of equal rates 
pi = . . . = pi\i, with different methods than presented here. 

CASE B: Bernoulli jumps with blocking. Particles are labelled from right to left, 
so YN{n) < Y]\f^i{n) < . . . < Yi{n). Between time n — 1 and n each particle attempts to 
move one step to the right, but it is constrained not to overtake the particle to its right. 
Particle i moves with probability pi. The particles are now updated from right to left, so it 
is the updated position of the particle to the right that acts as a block. 

The evolution is generated from a family [C{k, n);k £ {1, 2, . . . A'^}, n G N) of independent 
Bernoulli random variables satisfying P(^(A:,n) = +1) = 1 — P{^{k,n) = 0) = pk, via the 
recursions Yi(n) = Yi(n — 1) + ^(l,n), and for A; = 2, 3, . . . A^, 

n(n) =min(yfc(n- 1) +e(/c,n),yfe_i(n)). 

In particular, y is a Markov chain on W'^ . 

The process Y has been investigated by Rakos and Schiitz in the context of a frag- 
mentation model. On shifting the i-th particle i positions to the left, Y corresponds to the 
discrete-time totally asymmetric simple exclusion process (TASEP) with sequential updat- 
ing. Moreover, the process arises in a directed first-passage percolation model known as the 
Seppalainen model [18] with 'origin' (0, 0), where Y{n) corresponds to the vector of instants 
at which the sites {n,0), . . . , (n, N — 1) become wet. 

An important process arises if we scale Y as in the law of small numbers, i.e., by setting 
Pi = oii/M and considering Y{\_Mt\) for t G R+ as M — > oo. In the case of equal rates, the 
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resulting continuous-time Markov process describes, after a deterministic shift, the positions 
of particles in the (continuous-time) TASEP. This is the framework originally studied by 
Schiitz [T7], who derives the transition kernel of this process. It has recently been extended 
to particles hopping at different rates by Rakos and Schiitz [16]. 

The same continuous-time process is also of significant importance for series Jackson 
queueing networks as well as for a corner-growth model. In the queueing context, it repre- 
sents the cumulative number of departures from each of the queues; see [5]. In the corner- 
growth model, it represents the height of the first N columns. We refer to Konig's survey 
paper [11] for these and further connections, such as the relation between the Ath com- 
ponent of this process and directed last-passage percolation with exponentially distributed 
travel times. 

CASE C: Geometric jumps with blocking. Once again particles are labelled from 
right to left, so Y]\f(n) < Yi^-i{n) < . . . < Yi(n). Between time n — 1 and n each particle 
attempts to move a geometrically distributed number of steps to the right, starting with the 
leftmost particle. As in case B, a particle is blocked by the closest particle to its right if it 
tries to overtake another particle, but this time it is the old position of that particle that 
acts clS cl block. 

The evolution is generated from a family (^(A;, n);k£ {1, 2, . . . A}, n € N) of independent 
geometric random variables satisfying F{S^{k, n) = r) = (1 — Pk)Pk foi' ?^ = 0, 1, 2, . . ., via the 
recursions Yi{n) = Yi{n — 1) + ^(1, n), and for A; = 2, 3, . . . A^, 

n(n) =min(yfc(n- 1) +e(A:,n),yfc_i(n- 1)). 

As laid out by Draief et al. [6j, this recursion arises in the study of so-called tandem stores 
in series. The random vector Y represents the cumulative demand met at each of the stores 
when the first store is saturated. Alternatively, Y can be interpreted as the minimum-weight 
vector of certain weighted lattice paths. 

CASE D: Bernoulli jumps with pushing. Now it is natural to label particles left 
to right, so Yi{n) < Y2{n) < ... < Y]\f{n). We update from right to left, and preserve 
the ordering by pushing particles to the right. The evolution is generated from a family 
[^{k,n);k G {1, 2, . . . A^}, n G N) of independent Bernoulli random variables satisfying 
P{^{k, n) = +1) = 1 — P(^(A:, n) = 0) = pk, via the recursions Yi(n) = Yi{n — 1) + ^(1, n), 
and for A; = 2, 3, . . . A, 

Yfc(n) = max(yfc(n - 1) + i{k,n),Yk-i{n)) . 

The process Y is the discrete-time analogue of a particle system studied by Alimohammadi 
et al. [Ij, which has been studied further by Borodin and Ferrari 3J. It also plays an 
important role in the directed /ast-passage analogue of the Seppalainen model. 

We summarize the description of the four cases in the following table. 



Case 


A 


B 


C 


D 


Jump distribution 


geometric 


Bernoulli 


geometric 


Bernoulli 


Interaction 


pushing 


blocking 


blocking 


pushing 


Updating 


from left 


from right 


from left 


from right 



We need some well-known symmetric functions in order to present the Markov transition 
kernel of Y in each of the four cases. The rth complete homogeneous symmetric polynomials 
in the indeterminates ai, . . . otv is given by 

hj.{a) = a^^a2^---a^. 

fci>0,...,fe]v>0;fci+fc2H l-fc]v=r 

By convention /iq = 1 and = Q for r < 0. Now for 1 < i < j < A, let h^^\a) = hr{a^'^-^'^) 
where afe') is the A^- vector (0, . . . , 0, Oj+i, a^-^-2, . . . , aj, 0, . . . 0) obtained from a by setting 
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the first i weights, and tlie last N — j weights equal to 0. Equivalently it is the rth complete 
homogeneous symmetric polynomial in the indeterminates aj+i, . . . ,aj. We set hr''\a) = 
l(r = 0). 

We also need e^, the rth elementary symmetric function defined as 

kl<k2<---<kr 

In analogy with the complete homogeneous symmetric functions, we use the conventions 

ei^^\a) = l(r = 0) and bq = 1. We also set e^'^(a) = e,.(Q;(*-')) and = for r < 0. 
Given a function / on Z and a vector a € R^, we write 



(1) fL''\k) 

and 



Xr=o4'\a)f{k + i) iii<j, 



" ^ ' lEr=o#^(«)/(fc-^) if^<j, 

provided all scries converge absolutely. Our main theorem uses this notation for functions 
/ belonging to two different families, [wn',n € Z) and {yn'-,n G Z), for which the desired 
convergence holds. These families are defined through 

= > 0;n > 0) and = Q 1(0 < A; < n; n > 0). 

We write Wn^2 for fa^^ with / = Wm and define Wn^^^ ^n?2, and Vn^2 similarly. Moreover, we 
abbreviate the vector . . . iP^^l) by jo^^, and define the vector vr through tTj = Pi/ {1—pi). 

Theorem 1. T/ie transition kernel Qn of the process Y is given by the following expressions. 

CASE A: Geometric jumps with pushing. We have for y, y' G , 

N 

Qniy, y') = n [(1 - P^Tpi''"] {^-. ivi -yj + i-j)}. 

k=l 

CASE B: Bernoulli jumps with blocking. We have for y, y' G , 

N 

Qniv, y') = n [(1 - P^r^''"] det {uli^K^i -yj-i+j)}. 

k=l 

CASE C: Geometric jumps with blocking. We have for y,y' € , 

N 

Qn{y, y') = n [(1 - PkTpf''"] det {w^M -yj-i + j)]- 

k=l 

CASE D: Bernoulli jumps with pushing. We have for y, y' G , 

N 

Qn{y, y') = n [(1 - Pkr4''"] det {y'i -y, + i-j)}. 

k=l 

The remainder of this paper is devoted to a proof of this theorem. 
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3. The RSK correspondence and its variants 

In each of the four cases considered in the previous section, the Markov process of interest 
(Y{n);n G N) is constructed from a family (^^{k,n);k € {1, 2, . . . , A^}, n S N) of random 
innovations. In this section we wih construct a second Markov process Z from the same 
innovations data, using the RSK algorithm or one of its variants. In each case, we will 
be able to find the (Karlin-McGregor type) transition semigroup of Z and show that it is 
intertwined with the transition semigroup of Y. 

We give some definitions in order to describe the RSK type algorithms in the form we 
need. A partition A with k parts is an integer vector Ai, . . . , satisfying Ai > A2 > . . . > A^. 
Consider an array of strictly positive integers T = {Tij;! < i < k,l < j < Xi) oi shape 
A satisfying Tij < Tj+ij and Tij < Tjj+i (interpret Tij as infinity if it is undefined). We 
write T E if the integers in T do not exceed n and increase strictly down the columns, 
while A consists of at most parts. In the terminology of enumerative combinatorics, 
consists of semi-standard Young tableaux (SSYT) with at most N rows and content 
{1, . . . ,n}. Similarly, we write T G T^*" if the integers in T do not exceed n and increase 
strictly along the rows, while A consists of at most parts. We write sh (T) for the shape of 
T, which we consider to be an element of by padding the vector with zeros if necessary 
(N is fixed throughout). 

We study four different ways to associate a T^'^-valued process {V{n);n > O) to the 
data [^{k,n)]k G {!,..., A^},n G N), each corresponding to a different variant of RSK. 
More details on the different variants can for instance be found in [TJIH]. 

We begin by noting two methods for constructing a two-line array of the form 

f ai a2 as 
V h 62 63 

from the innovation data, where the are nondecreasing. Both methods have the property 
that the column appears ^{b, a) times in the array. The first method, lexicographic 
array- construction, requires that bi < ftj+i if Oj = Oj+i. The second method requires that 
bi > if ai = aj+i, and we therefore call it the anti-lexicographic array- construction. 
Next we describe two methods for constructing a sequence {V{n);n > O) of SSYT in 

T^^'^ from the given two-line array. Both constructions are inductive, and start with an 
empty SSYT 'P(O). Given V{n), the SSYT V{n -\- 1) is found by inserting the elements 
bi for which Oj = n + 1. If there are M such elements, the methods construct a sequence 
V^{n), . . . ,P*^(n) such that V^in) = V{n) and V{n + 1) = P*^(n). The SSYT V'+^{n) is 
constructed from V^{n) by inserting the next unused 6-element of the two-line array. The 
first method, row insertion, inserts an element b into V'^~^^{n) using the following rules: 

• If every entry in the first row of 'P^{n) is smaller than or equal to b, then b is appended 
to the end of the row. 

• Otherwise, b is used to replace the leftmost entry in the row which is strictly larger 
than b. 

The entry replaced is inserted in the same manner into the second row, and this process 
continues until an entry is either placed at the end of a row or it is placed in the first 
position of an empty row. The second method, column insertion, is a modification of the 
above procedure. Instead of inserting entries along the rows, the entries are now inserted 
down the columns. Now the following rules are followed: 

• If every entry in the column is strictly smaller than b, then b is appended to the end 
of the column. 

• Otherwise, b is used to replace the uppermost entry in the column which is greater 
than or equal to b. 

Four combinations of ^-values, array constructions, and insertion algorithms are of special 
interest, as there is a combinatorial correspondence underlying the above construction of 
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the process V. To explain this, let Q{n) be the unique array of integers for which the entries 
1, . . . ,m form an array with the same shape as V{m) for m < n. Depending on the RSK 
variant chosen, the pair (V{n), Q{n)) belongs either to 

(2) {{S, T) G T;^'^ X T^^ : sh {S) = sh (T)} 
or to 

(3) {(5, T) G T;^'^ X T^< : sh {S) = sh (T)}. 

The four combinatorial correspondences map the sequence (^(A;, m); k G {1, 2, . . . , A^}, m G 
{1, 2, . . . , n}) bijectively to ('P(n), Q{n)) for any n > 1. We present these four combinations 
in the next table, along with the name under which the resulting bijection (correspondence) 
is known. 



Correspondence 


RSK 


dual RSK 


Burge 


dual Burge 


Innovations data 


0,1,2,... 


0,1 


0,1,2,... 


0,1 


Range (7'(n),Q(n)) 


m 


(H 


© 


m 


Insertion algorithm 


row 


column 


column 


row 


Array construction 


lexicographic 


lexicographic 


anti- lexicographic 


anti-lexicographic 



It is our next aim to relate the four correspondences to the interacting-particle framework 
of Section [2j Further analysis is facilitated by a second process Z = [Z{n);n > O) arismg 
from the process V. We define this process by letting Z{n) be the shape of V{n), so Z 
takes values in W'^ . After giving the transition kernel of Z, we prove that the transition 
semigroup of Y is intertwined with the transition semigroup of Z. 

The semi-standard Young tableaux in xj^'^ play an important role in our analysis since 
the process V is xj^'^-valued, and we need some further definitions for such tableaux. For 
T G T^'^, we let ledge(T) be the A^-vector for which element i is the number of i's in row 
i (the terminology 'ledge' is motivated in the next section). This vector may contain zeros 
and we always have ledge(T) G W'^ . Similarly, we let redge(r) be the A^-vector for which 
the i-th entry is the number of elements in the first row that do not exceed i; we always 
have redge(T) G . For a vector a = (ai, a2, . . . , a^) of weights, we define the weight 
of T G T;^'^ by 

N 
i=l 

where ni[T) is the number of Vs in T. This is the usual definition for the weight of a SSYT. 

CASE A: Geometric jumps with pushing. Consider the RSK algorithm, i.e., lexico- 
graphic array construction and row insertion. The row-insertion algorithm shows that the 
vector- valued process (redge(P(n)); n > O) is exactly the same as the process (y(n);n > O). 
Note also that Y is Markov relative to the filtration of the Markov process V. 

The dynamics of the RSK algorithm show that Q{n) is a SSYT. Therefore, using the 
bijective property of RSK we obtain for S G T^'^, T G T^^, 

TV 

(4) V{V{n) = S, Q(n) = T) = J] [1 - p,,]"/ l(sh(5) = sh(r)). 

k=l 

Using the fact that Q(n) encodes the shapes (Z(m);m < n), as argued in [12, Sec. 3.2] 
we can use dH to compute the law of the shape process Z. Indeed, summing (jj]) over 
appropriate S G xj^'^, T G we obtain for m < n, 



k= 



N 

P{Z{n) = z{n), Z{m) = z(m\ . . . , Z(l) = z(l)) = J] [1 " VkT ^z(n){p)9:[j;'/,^^), 
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provided the left-hand side is nonzero. Here Sz{p) = YlxeT^''^- sh{T)=zP'^ ^ symmetric 
function (in p) known as a Schur polynomial. Also, dx/^ is the number of skew SSYT 
with shape X/fj, and entries from {1, . . . ,k} [20l Sec. 7.10]. The Jacobi-Trudi identity pOt 
Thm. 7.16.1] implies 

9\/f, = det {wk{Xi - j)} , 

and we thus find that Z is a Markov chain on with transition kernel given by 

Pn{z,z') = TT [1 -Pfc]" ^^T^det {Wn{z'i - Zj - i + j)} . 

k=l ^""^P' 

A similar reasoning, now summing @ over S E T^'^ with redge(5) = y € W'^ and 
sh(5') = Z{n), shows that 

P(y(n) = y\Z{m),m < n) = kp{Z{n),y), 

where 



- JV 



TeT^-'';sh{T)=z,redge{T)=y 



Note that P{Y(n) = y\Z{m),m < n) depends on (Z(m);m < n) only through Z{n). This 
yields for y G and m < n, 

^{Yin) =y\Z{k),k <m) = P [F{Y{n) = y\Z{k),k < n)\ Z{k),k < m] 

= p\^Kp{Z{n),y) Z{k),k <m 

= E Pn~rn{Z{m),z)kp{z,y). 

On the other hand, the left-hand side can be written as 

P{Y{n) =y\Z{k),k <m) = P [P{Y {n) = y\V{k),k < m)\ Z{k),k < m\ 

= P[Q n—m 

(y(m), y)\ Z{k), k <m] 

On combining these two displays, we deduce the intertwining relationship PnKp = KpQn, 
where the product PK of the kernel P with domain x and the kernel K with 
domain W'^ x is defined as PK(z,y) = "^^i^^tn P{z, z')K{z\y). The next section 
investigates the intertwining relationship in detail to find the kernel Qn- 

CASE B: Bernoulli jumps with blocking. Under lexicographic array construction 
and column insertion, (ledge(P(n)); n > O) is exactly the same as the process {Y{n)-n > 0). 
The bijective property of the dual RSK shows that for S G Tj^'^,T G T^*~, 

AT 

P{V{n) = S, Q{n) = T) = J] [1 - p^T vr^ l(sh(5) = sh(r)). 

fc=i 

As in [121 Sec. 3.3], on combining this with the dual Jacobi-Trudi identity, we find that Z 
is a Markov chain on with transition kernel 

P„(z, z') = n [1 - PkT ^ det {vn{z[ -z,-i + j)]. 



After setting 



^ ' TGT^.^;sh(T)=z,ledge(T)=y 



T 

a , 
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we may derive the intertwining PuKj^ = Kj^Qn along the hnes of case A. 

CASE C: Geometric jumps with blocking. The dynamics of the column insertion 
algorithm show that (ledge('P(n)); n > O) is exactly the same as the process (Y{n);n > O). 
Moreover, since the law of the process Z is invariant under the choice of the insertion 
algorithm, Z is Markovian with the same kernel as in case A. The intertwining is PnKp = 

CASE D: Bernoulli jumps with pushing. Now (redge('P(n)); n > O) is exactly the 
same as the process (y(n);n > O) . The process Z is Markovian with the same kernel as in 
case B, and we have the intertwining PuK^^ = Kj^Qn- 

4. Determinantal intertwining kernels 

It is the aim of this section to show how the kernel Qn of Y can be recovered from any 
of the intertwining identities PnKa = KaQn and PnKa = KaQm where Pn is the (known) 
Karlin-McGregor type kernel of the process Z. In fact, we prove the stronger assertion that 
the intertwining kernels Ka and are invertible. 

In what follows, it is convenient to embed T^'^ in a space parametrized by an array of 
variables x = (x^, . . . , x^) with x'^ = (x^, . . . , x^) G Z'^, such that the coordinates satisfy 
the inequalities 

Xk ^ — ^k-1 ^ ^k~2 ^■■■21X2^X1 ^ Xi 

for k = 2, . . . , N . The pattern x corresponding to a given SSYT in xj^'^ is found by letting 
3cj mark the position of the last entry in row i whose label does not exceed j. We write 
for the set of all x satisfying the above constraint, and say that any x € is a Gelfand- 
Tsetlin (GT) pattern. In contrast to the tableau setting, it is not required that x has non- 
negative entries. For x G K^, we set sh(x) = {x^ ,X2 , ■ ■ ■ ,x^), ledge(x) = {x\, . . . ,x^), 
and redge(x) = {x\, . . . ,x^); these definitions are consistent with those given for a SSYT 

in T^'^. For a vector a = (ai, 02, • • • j oin) of weights, we define the weight of a GT 
pattern x by 

N , 
a =«i ii^fc 

k=2 

in accordance with our previous definition for tableaux in T^^'^. 

Instead of studying and , it is equivalent but more convenient to work with modified 
versions which do not involve Schur functions but have a polynomial prefactor. It is natural 
to start with the kernel Ka, since this is a 'square' W'^ x matrix. We define, suppressing 
the dependence on a, 

(5) A(z,y)=a-^^ J] 

xeK^; sh(x)=^, ledge(x)=j/ 

We first show that A(z, y) can be written as a determinant. 

Proposition 2. For y,z £ , we have 

A{z,y) = det{hi^^^l_^^^{a)}. 

Proof. The proof is a variant of a well-known argument using non-intersecting lattice paths 
to derive the Jacobi-Trudi formulae, see for example ^Oj. Fix z and y belonging to W'^ 
with zn = UN- Each Gelfand-Tsetlin pattern x G having sh(x) = z and ledge(x) = y 
can be encoded as a non-intersecting (N — l)-tuple of paths (Pi, P2, . ■ ■ , Pn-i) on the edges 
of the square lattice with vertex set Z^. Paths always traverse edges in the direction of 
increasing co-ordinates: either 'upwards' or 'rightwards'. The path P^ begins at the vertex 
iVk — k,k + 1), ends at {z^ — k,N), and contains the 'horizontal' edges from (x^fT^ — k, r) to 
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(x^ — k, r), for r = k + l,k + 2,...,N. This correspondence between patterns and paths is 
a bijection and consequently we may write the sum defining A{z, y) as 

(Pl,...,PjV_l) 

where the sum is over all non-intersecting paths which have starting points and end points 
given in terms of y and z as above, and where the weight w{Pk) is defined to be 

N 

n ^r^^ 

r=k+l 



with e(r) denoting the number of horizontal edges at height r contained in the path P^. 
By the Gessel-Viennot formula (see Theorem 2.7.1 of [19]), this sum for A{z^y) is equal to 
det(M) where M is an [N — 1) x [N — 1) matrix with (i, j)th entry given by w{P) where 
P runs through all paths connecting — j, j + 1) to {zi — i, N). This latter quantity is easily 

found to equal h^^^'^y._i_^_j{a). The proposition follows on noting that h^J^_f^J^_-j^]^{a) = 
for i = 1, 2, 3, . . . , — 1 and that h['^^y^{a) = l{zj\f = yN)- D 

We will use the following identity that is easily checked by means of generating functions, 

UN), X ^ J /in*-* (a) if i < ^> 



(6) 



Recall also the Cauchy-Binet formula 

(7) ^ det{<l)i{zj - j)}dei{iij{zi- i)] = dei}/Y^(j)i{z)^j( 

Proposition 3. The x matrix A is invertible. Its inverse is given by the x 
matrix H defined as 

U{y,z) = det{{-l)y^-^^-'+^ef_l^_,^^{a)}. 

Proof. We first show that 11 is a left inverse of A. From the Cauchy-Binet formula we deduce 
that 

(8) nA(y,y') = det{m^^,_^^^}, 



where 




(ii) 

Observe that mi- is zero if either i < j and r > j — i ot ii i > j and r < 0. In particular 
the product of the diagonal elements appearing in the determinant on the righthandside of 
([8]) is l{y[ = yi). We need to show that this is the only contribution to the determinant. 
Suppose that vr is a permutation of {1,2, . . . , A^} other than the identity, and consider the 
product of the (i,7r(«))th entries. It easy to see that there must be integers ii and 12 such 
that ii < i2, TT{ii) > ii, 7r(i2) < ^2 and 7r(zi) > 7r(z2). Now the (ii, 7r(zi))th entry is zero 
if > y^(j^), and the (i2, vr(i2))th entry is zero if yi^ < y'^f^^^y But since i/i, > yi^ and 
y^(j^) < 2/^r(i2)' least one of these previous inequalities holds, and hence the product is 
zero. 

The key ingredient for showing that II is a right inverse of A is the trivial identity 
U{y,z')= Y: (-l/^4°^)..•(-l)^-4;^^)det{l(y.-^ = 4-J+^.)}. 

£l,.../jvGZ 



10 A. B. DIEKER AND J. WARREN 

Observe that the sum over the i is finite since the summand is zero unless < ii < N — i 
for all i. We therefore have for any function / on and < £i < N — i, 

f{yi -l,...,yN-N) det{l{yi - i = - j + 
= J2 sgn(^)/(4(i) - a{N) + £n) 

N 

X XI n ^(y™ -m = z'^^,^^ - a{m) + C) 

ygVK^ m=l 

= f{z[-l+ei,...,z'^-N + £N), 
where Sn is the set of permutations on {1, . . . , A^}. After absorbing the sum over the i in 
the determinant, we obtain AIl{z,z') = det{m^^''^^,_._^^} = l{z = z'). □ 

By virtue of this proposition, the intertwining AQ = PA yields Q = IIPA, and a straight- 
forward computation using the Cauchy-Binet formula and ([6]) implies the following corollary. 
Recall the definition of fa^^ in ([1]). 

Corollary 4. Suppose that P is a x W'^ matrix of the form 

P{z,z') = det {/(z- - Zj - i+j)} 

for some function f on Z. Suppose that Q is another x matrix and that the 

intertwining relation PA = AQ holds. Then we have for y,y' € W'^ , 

Q{y,y') = det{f<'^\y',-y,-i + j)]. 

We will also consider a second intertwining kernel that arises by replacing the left edge 
of the pattern with the right edge in ([5]): for z G , y G , we set 

A(z,y)=ar^i X 

xeK-*^; sh(x)=2, redge(x)=j/ 

For X € K^, define x G by = —x'l_-j^-^. The correspondence x i— >• x is bijective and it 

is easily verified that = where (3 = = (aj~^, . . . , a]^^)- Using this correspondence 
and our results for A we obtain the following. 

Proposition 5. For z G and y G , we have 

Hz,y) = det{h^^X-...^a^"^}- 
Moreover, A is invertihle with inverse 11 given by, for y G and z G W'^ , 

U{y,z) = det{(-l)^---^-^-^+^et^_)^^^_^^„,^^.(a-i)}. 

The analogue of Corollary [J] follows immediately from this proposition. 

Corollary 6. Suppose that P is a x VF^ matrix of the form 

P{z,z') = det {/(z- - Zj - i+j)} 

for some function / on Z. Suppose that Q is a x W'^ matrix and that the intertwining 
relation PA = AQ holds. Then we have for y, y' G , 

Q{y,y') = dei[f^%[-y,+i-j)]. 
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